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Abstract 

First-principles calculation predict that olivine Li4MnFeCoNiP40i6 lias ferro- 
toroidic characteristic and ferrimagnetic configuration with magnetic moment of 
per formula unit. The ferrotoroidicity of this material makes it a poten- 
tial candidate for magnetoelectric materials . Based on the orbital-resolved density 
of states for the transtion-metal ions in Li4MnFeCoNiP40i6, the spin configura- 
tion for Mn2+,Fe3+,Co2+, and Ni2+ is t%{])el{^) ti,(i)e2(i) 4,(i)ti,(T)ei (i)e2(T) 
and i2g(T)^2g(i)6g(T)eg(i) , respectively. Density functional theory plus U (DFT-I-U) 
shows a indirect band gap of 1.25 eV in this predicted material, which is not simply 
related to the electronic conductivity in terms of being used as cathode material in 
rechargeable Li-ion batteries. 
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1 Introduction 



Ferrotoroidic(FTO) domains liave been observed in olivine LiCoP04 recently 
using second harmonic generation (SHG), which are independent of the anti- 
ferromagnetic (AFM) domains [1]. The ordered arrangement of magnetic vor- 
tices is found to be the fourth form of ferroic order, which is currently debated 
whether to be included or not [2|3Pf5|6] . The well-known three forms of ferroic 
are ferromagnetism, ferroelectricity, and ferroelasticity. This fourth form of 
ferroic is named as ferrotoroidicity corresponding to the former three forms 
,which have been established well and investigated extensively. A ferrotoroidic 
moment is generated by a vortex of magnetic moment, which is asymmet- 
ric under the reversal of time and space. The time- and space- asymmetric 
property relate ferrotoroidic materials deeply to the multiferroic materials, 
in which ferromagnetic property violates time reversal symmetry and ferro- 
electric property violates spatial inversion symmetry. In multiferroic materials, 
magnetization can be rotated by an external electric field, and polarization can 
be reversed by a reversal magnetic field, which is called the magnetoelectric 
effect. Therefore, FTO materials can be used as magnetoelectric materials, in 
which a magnetization is induced by an electric field, and vice versa, and have 
potential applications in information storage, the emerging field of spintronic 
, and sensors. 

In olivine LiCoP04, the Co^+ ions are located at positions (l/4+e,l/4,-(5), 
where e and 6 represent the small displacements by the mmm symmetry. The 
four nearest Co^^ are arranged antiferromagnetically along the y axis and 
have two different toroidic moments due to the rotation of spins by </> = 4.6° 
away from y axis[7j.The two toroidic moments are opposite in direction and 
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not equal in value, leading to the ferrotoroidic porperty. Magnetism can be 
improved by substitution with the transition-metal ions by another one such 
as in Bi2FeCr06[8j and Bi2FeTi06^. In this paper, we put the high-spin 
Mn^+ and Fe^"*" in the sites with larger toroidic radius while low-spin Co^+ 
and Ni^"*" in the sites with smaller toroidic radius, which would bring about 
a large resulting toroidic moments and further an excellent ferrotoroidicity. 
This predicted olivine structural material, Li4MnFeCoNiP40i6, is expected to 
have significant magnetoelectric effect as multiferroic materials due to ferro- 
toroidicity, and also have potential to be used as cathode materials for Li-ion 
rechargeable batteries. In order to gain much insight in this predicted ma- 
terial, we have investigated the electronic and ferrotoroidic properties using 
first-principles calculations. The remainder of this article is organized as fol- 
lows: in section 2, we present the computational details of our calculations. In 
section 3, we report the calculated results and discussion. In section 4, the 
conclusions based on our calculation are given. 

2 Computational details 

Calculations in this work have been done using the Quantum-ESPRESSO 
package [TU]. Our density functional theory (DFT) computations follow the 
method previously reported in Ref. ^.We used our self-interaction-corrected 
ultrasoft pseudopotential implementation with the Perdew Burke Ernzerhof 
(PBE) exchange correlation functional, as the common local density approxi- 
mation(LDA) fails to obtain a band gap in the transition-metal oxides. Plane- 
wave basis set with kinetic energy cut-off of 38 Ry was employed. Li 2s, Mn 
3s, 3p and 3d electrons, Fe 3s, 3p and 3d electrons, Co 3s, 3p and 3d elec- 



3 



trons, Ni 3s, 3p and 3d electrons, P 3s and 3p electrons, and O 2s and O 2p 
electrons have been treated as valence states. We used up to 6 x 7 x 8 grids of 
special k points in total energy calculation, and 8 x 9 x 10 grids of special k 
points were used for density of states(DOS) calculation. The lattice constants 
and the atomic positions were taken from Ref . [H] and [12] . Then, fully struc- 
tural relaxations were performed by minimization of the Hellman-Feynman 
forces within a convergency threshold of 10"'^ Ry/Bohr. The transition-metal 
ions were arranged differently to find the favorable structure with minimal 
energy. Consequently, the stable structural configuration is obtained and con- 
structed by XCrySDen[13j as in Fig.l, in which the unit cell doubled along 
the b axis. Collinear spin configurations have been used to construct the dif- 
ferent AFM order. The favorable configuration is obtained as shown in Fig. 

2 accordingly. Local density approximation(LDA) or generalized gradient ap- 
proximation(GGA) within DFT always predicts the conduction properties in 
transition-metal oxides, and tends to over-delocalize electrons when the ki- 
netic energy of electrons is not large enough to overcome the on-site repulsion. 
Hence, a Hubbard like term U is introduced and taken as in Ref.[T5 to treat 
the transiton-metal ions, which is named as DFT+U. 

3 Results and discussion 

It is known that LiMP04(M=Mn, Fe, Co, and Ni)has an orthorhombic cell, 
which contains four formula units, 28 atoms: 4 lithium, 4 M, 4 phosphor, 
and 16 oxygen atoms. Phosphor ions are tetrahedrally coordinated to oxygen, 
where M and lithium ions occupy the centers of distorted oxygen octahedral 
with space group Pnma as shown in Fig.l. The lattice parameters and atomic 
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positions are from Ref.[TT] and [12], and we put the Mn, Fe, Co, and Ni ions on 
the 4M sites, respectively. The stable structure under the condition of minimal 
total energy is shown in Fig. 1. The high-spin Mn^"*" and Fe^^ occupy the sites 
having large toroidic radius, which tends to generate a net toroidic moment 
as we expected. In Fig. 2, atoms in one unit cell are projected on the XY 
plane to make clear the relative positions of these four type of ions. Next, the 
lattice parameters and atomic positions are fully relaxed. The relaxed lattice 
parameters and Wyckoff positions are shown in Table 1. The lattice parameters 
are about the averaged values of LiMP04(M=Mn, Fe, Co, and Ni) due to the 
mixing of different transition-metal ions. 



Let's discuss the magnetic ordering of these four type of ions in the system. The 
spin configurations are constructed as no spin, ferromagnetic, and four anti- 
ferromagnetic configurations(denoted as NM, FM, AFMl, AFM2, AFM3 and 
AFM4, respectively.). The AFM configuration are shown in Table 2 in detail, 
and the four atomic sites to different AFM configuration are Co,Mn,Ni,and 
Fe as shown in Fig. 2 anticlockwise, and the y axis is taken as the positive di- 
rection. The differences in the total energy in the NM, AFMl, AFM2, AFM3, 
and AFM4 configuration relative to its total energy in the FM configuration 
are shown in Table 2. One can see that the NM configuration has the highest 
total energy and becomes the most unstable configuration. The AFM2 is the 
optimum magnetic configuration shown in Fig. 2. A resulting magnetic mo- 
ment of 1.56 /is is obtained from calculation, which is related deeply with 
the real spin configuration of the four transition-metal elements, which will be 
discussed in detail in the following sections. 
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The spin part of the toroidic moment is expressed as 

Toc^r„xs„, (1) 

n 

where r„ is the radius vector and s„ is the spin of the n magnetic ions, and 
the center of the unit cell is taken as the origin. As we know, the Mn and Fe 
ions are expected to be in the higher spin states as compared with Co and Ni 
ions. Therefore, as long as spins rotate away from y axis, from the favorable 
structure shown in Fig. 2, we get 

^Mn,Fe,n X SMn,Fe,n > ^Co,Ni,n X ^Co,Ni,n- (2) 

Hence, a large resulting toroidic moment can be obtained in this predicted 
system as compared with the case in LiCoP04. From the calculation above, 
we consider the ferrotoroidic and ferrimagnetic properties can be achieved 
simultaneously. 

The DFT+U method combines the high efficiency of LDA/GGA, and an ex- 
plicit treatment of correlation with a Hubbard-type model for a subset of 
states in the system. The total energy is expressed as 

Elda+uHt)] = EldaHt)] + EnubiW^}] - E^clW'}], (3) 

where n(r) denotes the charge density and n^^ is the transition-metal on-site 
occupation matrix. From the expansion in spherical harmonics, in first order 
approximation, we have 

Elda+u = ELDA[n{r)] + Eu[{nZrA] = Elda[p] + ^ E Tr[n'\l - n'')] (4) 

^ i,s 

In present work, U is 3.92, 3.71, 5.05, and 5.26 for Mn, Fe, Co, and NifH], 
respectively. 
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The total density of states(DOS) using DFT and DFT+U for Li4MnFeCoNiP40 
are shown in Fig. 3, and the DOS of LiC0P04 is also reported for compari- 
son. From the DOS of LiC0P04 and Li4MnFeCoNiP40i6, it is clear that the 
manifold states near the Fermi level in Li4MnFeCoNiP40i6 are attributed to 
the other transition-metal elements(Mn, Fe, and Ni) except Co. It can be seen 
that DFT predicts a half-metal behavior, in which there is only a complete 
spin polarization of electrons at the Fermi level. The strong localization of 
3(i-electrons of these transition-metal leads to the states passing through the 
Fermi enrgy, and further to the half-metal characteristic. The strong correla- 
tion of these 3d electrons are well described by employing the DFT+U, and a 
indirect band gap of ~1.25 eV is found in the predicted material. 

In Fig. 4 we present the DFT+U band structure of Li4MnFeCoNiP40i6for ma- 
jority spin. It can be seen that it is a semiconductor with an indirect band gap 
~ 1.25 eV. The valence band top and conduction band bottom lie, respectively, 
at S and F points of the first Brillouin zone. This suggests that the electronic 
conduction performance of this new material is not very good. A large band 
gap will lead to a very small number of intrinsically generated electrons or 
holes when it is used as the cathode material for Li-ion rechargeable batteries. 
Therefore, the band gap will not play any significant role in setting the con- 
centration of conduction electrons or holes. It is likely that the key electrons 
involved in transport in the predicted material are not delocalized electrons, 
but localized small polarons[T5], and polaron mobility is determined by the 
hopping rate of the polarons. On the other hand, an excellent magnetoelectric 
effect can be achieved by the insulating property , in which the current leakage 
is well inhibited. Moreover, from the band dispersion, the valence bands are 
very dense with little dispersion, while the conduction bands have distinctive 
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dispersions, and these dispersive conduction bands are well above the Fermi 
energy. 

In order to gain insight in the spin configuration and the oxidation states 
of the transition-metal elements in the new predicted system, we report the 
orbital-resolved DOS for Mn, Fc, Co, and Ni in Fig. 5, Fig. 6, Fig. 7, and 
Fig. 8, respectively. From Fig. 5 , the majority-spin(up-spin) Mn-3d states in 
dxy, dyz, dxz, and dr^2_y2 orbitals are all occupied, and minority-spin (down- 
spin) in corresponding orbitals are all empty. This result is the formal high- 
spin(S=5/2) state of Mn^"*" with ^2g(t)Cg(t) configuration. Moreover, the majority- 
spin in (1^2 orbital is pushed away from the Fermi level as compared with the 
other orbitals, which is partly caused by the specific crystal field. It is worth 
pointing that all minority-spin are well above the Fermi level. 

From Fig. 6, the majority-spin Fe-3rf states in dxyidyz-,dz2,dxz,Ska<\ dx2-y2 or- 
bitals are all occupied, and minority-spin in corresponding orbitals arc all 
empty. The partial DOS of Fe-3(i indicates the high-spin(S=5/2) state of Fe^"*" 
with t\g{[)e^g{[) configuration. The majority- and minority-spin are all pushed 
rather away from the Fermi level due to applying the DFT-I-U, which improves 
the correlation of 3(i-electrons. The minority-spin in d^^ orbital is far away as 
in Mn-3(i states. 

From Fig. 7, majority-spin for Co-3d states in all five orbitals are filled, while 
minority-spin in dxy, dxz, and dx2-y2 are partially occupied, and in dz^ are fully 
filled, indicating the complicated spin configuration. Based on the partially 
occupied three orbitals and fully occupied one orbitals in minority-spin states, 
we regard the Co as the high-spin(S=3/2) state of 00^+ with tlg{i)tlg{])e]{i 
)eg(|) configuration. 
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Prom Fig. 8, It can be seen the spin configuration for Ni is rather comphcated 
as compared with the other transition-metal elements in the system. Majority- 
and minority-spin are all filled in Ni-3d states. Majority-spin in dr^i-y2 and 
dxz are near the Fermi level, which is partly caused by the employment of 
DFT+U, and some of these majority-spin would lie above the Fermi level 
within DFT. Based on the occupation of all spin states, we consider that Ni 
is in low-spin(S=2/2) state of Ni^+ with configuration. Interestingly, the 
three majority-spin electrons have possibility to occupy all five orbitals, and 
this point need to be studied further. This feature would decrease the magnetic 
moment of Ni^"^, which is in good agreement with the net magnetic moment of 
1.56 /is calculated. Based on the assumption that the majority-spin near Fermi 
level in d^p^y-i and d^z is derived from the conduction bands within DFT-I-U, 
we consider that the spin configuration for Ni^+ is ^2g(T)^2(/(J')^g(T)6g(i) 

Based on the structural analysis and spin configuration, the toroidic moment 
caused by Mn and Fe is greater than that caused by Co and Ni, which leads to a 
net toroidic moment and further to ferrotoroidic characteristic, accompanying 
with a ferrimagnetic configuration. This material can be used as multiferroic 
material due to its magnetoelectric coupling effect, and possible candidate for 
cathode material in rechargeable Li-ion batteries. For the latter, the electronic 
conductivity is not only simply related to the indirect band gap, but also 
related to localized small polarons and polaron conduction. 

4 Conclusions 

Motivated by the ferrotoroidic domains observed in LiCoP04, we predicted 

a new material by means of partial substitution of Co by Mn, Fe, and Ni 
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ions, in which the high-spin Mn and Fe ions are placed on the sites with large 
toroidic radius. Therefore, this material is expected to have excellent ferro- 
toroidictiy, and it can be used as multiferroic material in terms of its mag- 
netoelectric effect due to the ferrotoroidicity. The new system also possesses 
a ferrimagnetic configuration with magnetic moment of 1.56 fiB- DFT+U re- 
veal that the indirect band gap is ~1.25 eV, the electronic conductivity in 
Li4MnFeCoNiP40i6 is not simply related to the band gap. The valence state 
is found to be Mn^+,Fe^+,Co^+, and Ni^+ for Mn, Fe, co, and Ni ions, respec- 
tively. 
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Table 1 

Relaxed Itillicx^ parameters and wyckoIT |")()si(i()ns. 



a(A) b (A) c (A) V{Af Wyckoff positions 

10.3324 6.0105 4.6922 291.3992 Li 0.0000 0.0000 0.0000 

Mn 0.7787 0.2500 0.5210 
Fe 0.2213 0.7500 0.4790 
Co 0.2787 0.2500 0.9790 
Ni 0.7213 0.7500 0.0210 
P 0.0964 0.2500 0.4218 

01 0.0997 0.2500 0.7476 

02 0.4546 0.2500 0.2019 

03 0.1704 0.0415 0.2804 

Table 2 

AFM magnetic configuration and calculated differences in total energy {AE) in the 
NM, AFMl, AFM2, AFM3, and AFM4 configuration relative to its total energy in 
the FM configuration in units of eV/formular unit. The four atomic sites to different 
AFM configuration are Co, Mn, Ni, and Fe anticlockwise. The y axis is set as the 
positive direction. 

AFMl AFM2 AFM3 AFM4 NM 
magnetic configuration HH H]] 

AE -0.02 -0.08 -0.05 -0.07 1.35 
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Figure captions: 

Fig.l Crystal structure of Li4MnFeCoNiP40i6 showing two unit cells 
constructed by XCrySDenfRef . [T3] ) . 

Fig. 2 The favorable antiferromagnetic configuration projected onto the XY 
plane. The length of the arrow indicates the strength of magnetism 

Fig. 3 Total DOS of Li4MnFeCoNiP40i6 within (a)DFT+U ,and(b) DFT. 
(c) show the DOS of LiCoP04 within DFT+U. 

Fig. 4 Electronic band structure of Li4MnFeCoNiP40i6 for the majority-spin 
states. The doted line indicates the Fermi energy. 

Fig. 5 Orbital-resolved DOS for Mn in Li4MnFeCoNiP40i6. (a), (b), (c), (d), 
and (e) show the DOS for d^y, dyz, dz^, dxz, and dx2_y2 orbitals respectively. 
Majority-spin states are shown in the upper portions and minority-spin 
states in the lower portions in all panels. The Fermi level is set to zero. 

Fig. 6 Orbital-resolved DOS for Fe in Li4MnFeCoNiP40i6. (a), (b), (c), (d), 
and (e) show the DOS for d^y, dyz, dz^, d^^z, and dx2_y2 orbitals respectively. 
Majority-spin states are shown in the upper portions and minority-spin 
states in the lower portions in all panels. The Fermi level is set to zero. 

Fig. 7 Orbital-resolved DOS for Co in Li4MnFeCoNiP40i6. (a), (b), (c), (d), 
and (e) show the DOS for d^y, dyz, dz2, dxz, and dx2-y2 orbitals respectively. 
Majority-spin states are shown in the upper portions and minority-spin 
states in the lower portions in all panels. The Fermi level is set to zero. 

Fig. 8 Orbital-resolved DOS for Ni in Li4MnFeCoNiP40i6. (a), (b), (c), (d), 
and (e) show the DOS for d^y, dyz, dz:2, d^z, and dx2_y2 orbitals respectively. 



13 



Majority-spin states are shown in the upper portions and minority-spin 
states in the lower portions in all panels. The Fermi level is set to zero. 
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Fig. 1. Crystal structure of Li4MnFeCoNiP40i6 showing two unit cells constructed 
by XCrySDen(Ref.[I3]). 




Fig. 2. The favorable antiferromagnetic configuration projected onto the XY 
plane. The length of the arrow indicates the strength of magnetism. 
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Fig. 3. Total DOS of Li4MnFeCoNiP40i6 within (a)DFT+U , and(b) DFT. (c) show 
the DOS of LiCoP04 within DFT+U. 




Fig. 4. Electronic band structure of Li4MnFeCoNiP40i6 for the majority spin states. 
The doted line indicates the Fermi energy. 
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Fig. 5. Orbital-resolved DOS for Mn in Li4MnFeCoNiP40i6. (a), (b), (c), (d), and (e) 
show the DOS for dxy,dyz,dz2,dxz, and dy.2_y2 orbitals respectively. Majority-spin 
states are shown in the upper portions and minority-spin states in the lower portions 
in all panels. The Fermi level is set to zero. 
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Fig. 6. Orbital-resolved DOS for Fe in Li4MnFcCoNiP40i6. (a), (b), (c), (d), and (e) 
show the DOS for dxy,dyz,d^2,dxz, and dj.2_y2 orbitals respectively. Majority-spin 
states are shown in the upper portions and minority-spin states in the lower portions 
in all panels. The Fermi level is set to zero. 
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Fig. 7. Orbital-resolved DOS for Co in Li4MnFeCoNiP40i6. (a), (b), (c), (d), and (e) 
show the DOS for dxy,dyz,d^2,dxz, and d^.'^-y^ orbitals respectively. Majority-spin 
states are shown in the upper portions and minority-spin states in the lower portions 
in all panels. The Fermi level is set to zero. 
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Fig. 8. Orbital-resolved DOS for Ni in Li4MnFeCoNiP40i6. (a), (b), (c), (d), and (e) 
show the DOS for dxy,dyz,d^2,dxz, and d^^-yi orbitals respectively. Majority-spin 
states are shown in the upper portions and minority-spin states in the lower portions 
in all panels. The Fermi level is set to zero. 
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